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ABSTRACT 

We fit the data for the binary-lens microlensing event MACHO 
98-SMC-l from 5 different microlensing collaborations and find two 
distinct solutions characterized by binary separation d and mass ratio q: 
(d,q) = (0.54,0.50) and (d, q) = (3.65,0.36), where d is in units of the 
Einstein radius. However, the relative proper motion of the lens is very 
similar in the two solutions, 1.30 kms -1 kpc" 1 and 1.48 kms -1 kpc -1 , 
thus confirming that the lens is in the Small Magellanic Cloud. The close 
binary can be either rotating or approximately static but the wide binary 
must be rotating at close its maximum allowed rate to be consistent 
with all the data. We measure limb-darkening coefficients for five bands 
ranging from / to V. As expected, these progressively decrease with 
rising wavelength. This is the first measurement of limb darkening for a 
metal-poor A star. 

Subject headings: astrometry, gravitational lensing, dark matter 
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1. Introduction 

The binary- lens microlensing event MACHO 98-SMC-l, found by the MACHO 
collaboration (Alcock et al. 1999a) in observations toward the Small Magellanic 
Cloud (SMC), was observed by five different groups. Each group attempted to 
measure or constrain the relative lens-source proper motion by fitting the observed 
light curve to a binary-lens model, and each concluded that the proper motion is 
consistent with the lens being in the SMC rather than the Galactic halo (Afonso et al. 
1998 [EROS]; Albrow et al. 1999a [PLANET]; Alcock et al. 1999a [MACHO/GMAN]; 
Udalski et al. 1998b [OGLE]; Rhie et al. 1999 [MPS]). Despite this unanimous 
opinion, there are several reasons for taking a closer look at this event. 

First, Albrow et al. (1999c) have subsequently developed a more robust method 
for finding binary- lens solutions that are consistent with a given light-curve data set. 
They found a broad set of degeneracies for the fit to MACHO 98-SMC-l based on 
PLANET data and showed that based on these data alone, the proper motion could 
not be constrained to better than a factor of four. Such degeneracies are likely to 
be endemic to binary-lens fitting (Dominik 1999a). Albrow et al. (1999c) showed 
that by including additional data it would be possible to remove at least some of 
these degeneracies, but they argued that one major ambiguity (between "wide" and 
"close" binaries) might be difficult to resolve. Hence, it is important to determine 
whether these degeneracies can in fact be resolved by combining all available data. 
In particular, the OGLE and MPS data together constrain the first caustic crossing, 
the MACHO data constrain the baseline, the PLANET data give excellent coverage 
of the main part of the second caustic crossing, and the EROS data give similar 
coverage of the end of the second caustic crossing. 

Second, caustic-crossing binary events allow one in principle to measure the 
limb-darkening profile of the source star, as Albrow et al. (1999b) have done for a 
K giant using MACHO 97-BLG-28. The source in MACHO 98-SMC-l is an A star 
as determined both from its colors (Albrow et al. 1999a; Rhie et al. 1999) and its 
spectrum (Albrow et al. 1999a). Since it is in the SMC, it is almost certainly metal 
poor. If the limb darkening were measured, it would be the first such measurement 
for a metal-poor A star. The caustic crossing occurred while the source was visible 
from South Africa. The PLANET SAAO data have excellent coverage of the main 
part of the caustic crossing, and it therefore might be possible to measure the limb 
darkening from these data alone using a variant of the method of Albrow et al. 
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(1999c) which is described more fully in § 3. However, it is not clear to what degree 
the degeneracies in the overall fit would compromise such a determination. By 
fitting all the data these degeneracies could be partially or totally removed. The 
EROS data from Chile provide excellent coverage of the end of the crossing and the 
MACHO / GMAN data also cover the end of the crossing. Neither of these data sets 
can determine the limb darkening without fixing the characteristics of the caustic 
crossing, which in turn requires the PLANET data. 

Third, the light-curve coverage of MACHO 98-SMC-l is one of the best of any 
binary-lens event yet observed. It is therefore an excellent laboratory in which to 
search for additional, unanticipated anomalies that may be present in microlensing 
events but have not yet been noticed. 

In § 2 we briefly review the data that are available for this event. In § 3 we 
summarize and extend the Albrow et al. (1999c) method for finding binary-lens 
solutions. In § 4 we present our results for static binaries, including measurement of 
the limb-darkening coefficients. In § 5, we derive the proper motion which in turn 
determines the projected separation of the binary. We use the latter quantity to 
constrain the period of binary orbit. We consider rotating binary models that satisfy 
this constraint in § 6. Finally, in § 7, we study the relationship of the solutions 
derived here to those reported in earlier investigations that were based on subsets 
of our combined data set. We show that all the close-binary solutions are in fact 
different positions within one broad minimum in x 2 - The wide-binary solutions of 
Albrow (1999c) represent another broad minimum. We also resolve some puzzling 
discrepancies between different solutions. 



2. Data 

We describe all dates using HJD' =HJD-2450000.0 where HJD is the 
Heliocentric Julian Date. The reported times are the midpoints of the exposures. 
We combine a total of 14 light curves obtained at 8 different telescopes. These were 
reduced using various photometry packages as described below. In all cases, points 
that failed internal tests of these packages were eliminated prior to beginning the 
fitting process. 

The first two light curves were taken in the (broad non-standard) MACHO R 
and MACHO B filters on the 1.3 m telescope at the Mount Stromlo & Siding Springs 



- 6- 



Observatory (MSSSO) near Canberra, Australia. These contain 727 and and 735 
points respectively, beginning about 5 years before the event and ending 81 days after 
the second caustic crossing at HJD' ~ 982.6. The Mt. Stromlo 1.3 m is normally 
used to search for microlensing events and hence typically takes one exposure per 
clear night. However, because of the importance of this event, a total of 23 exposures 
were obtained during the nights just before and after the second caustic crossing. 
The next two curves were taken in the (standard Johnson/Cousins) R and B filters 
on the 0.9 m telescope at the Cerro Tololo Interamerican Observatory (CTIO) 
near La Serena, Chile. These contain 83 and 22 points respectively, beginning 7 
days before the first caustic crossing at HJD' ~ 970.5 and ending 6 days after the 
second crossing. MACHO has 1 hour per night on this telescope which was mostly 
dedicated to MACHO 98-SMC-l during this period of observation. All four of these 
light curves were reported by Alcock et al. (1999a). However, all the points after the 
first caustic crossing of the MSSSO data were re-reduced using image subtraction 
(Tomaney 1998) which we determined has somewhat smaller errors than the original 
SoDoPhot reductions. These include respectively 16 and 21 late-time points (more 
than 24 days after the second caustic crossing) that were not previously reported by 
Alcock et al. (1999a). 

Next is the OGLE (standard Cousins) I band curve from the 1.3 m Warsaw 
telescope at Las Campanas, Chile. The images are from OGLE's routine monitoring 
of the SMC to search for microlensing events, and OGLE made no special effort 
to observe this event. Rather, they analyzed their images after the event was over 
and found 7 measurements during the interval from 4 days before the first caustic 
crossing to 4 days after the second. As discussed by Udalski et al. (1998b), the 
primary interest in this relatively small data set comes from the second data point 
on HJD' = 970.9037 which is highly magnified (A ~ 29) and therefore comes either 
just before or just after the first caustic crossing. 

The sixth and seventh curves were taken in (broad non-standard) EROS R and 
EROS B filters on the 1.0 m Marly telescope at the European Southern Observatory 
at La Silla, Chile. Normally, this telescope is operated in survey mode to search for 
microlensing events. However, it was down for maintenance during most of the time 
that MACHO 98-SMC-l was inside the caustic, and recommenced operations only 
on the night of the second caustic crossing. In light of the importance of MACHO 
98-SMC-l, the telescope was entirely dedicated to observing this event during this 
night and made several observations per night for the next 15 nights, whereupon 



-7- 



it resumed normal operations. The observations of the first night were previously 
reported by Afonso et al. (1998). The rest of the observations are reported here for 
the first time. All the observations have been re-reduced using an improved version 
(Alard 1999) of the algorithm used by Afonso et al. (1998). There were a total of 
111 observations in R and 131 in B. These include about 8 points in each band from 
two years before the event and about another 8 from the year after the event when 
the source is approaching baseline. 

The eighth curve is based on the (standard Cousins) R band observations 
taken by the MPS collaboration using the 1.9 m telescope at MSSSO. A total of 
34 observations were taken from just before the first caustic crossing until 4 days 
after the second. In addition there is one late-time baseline measurement taken 67 
days after the caustic crossing. These observations were earlier reported by Rhie et 
al. (1999). Of particular note is the first observation on HJD' = 970.0485 which is 
clearly before the caustic crossing. This data point, combined with the OGLE data 
point 0.9 days later, strongly constrain the time of the first crossing. 

The ninth through thirteenth curves are based on standard Cousins / band 
observations taken by the PLANET collaboration using the SAAO 1 m telescope at 
Sutherland, South Africa, the Yale-CTIO 1 m, the CTIO 0.9 m, and the Canopus 
1 m near Hobart, Tasmania. The SAAO data are divided into two groups because 
of a change in the CCD detector at HJD' = 980.0, two days before the caustic 
crossing. The five PLANET data sets comprise respectively 13, 175, 32, 13, and 1 
observations. The CTIO 0.9 m data cover only the interior of the caustic, beginning 
3 days after the first crossing and ending 4 days before the start of the second. 
The Canopus data contain only 1 point about 1 day before the second crossing. 
The Yale-CTIO 1 m data begin 4 days before the second crossing and end 14 days 
after it. The SAAO 1 m data begin 5 days after the first caustic and end 44 days 
after the second. Of all the observations, only the SAAO data cover the peak of 
the second caustic crossing. Moreover, they do so quite densely. Most of these 
data were previously reported by Albrow et al. (1999a, 1999c). However, all of the 
SAAO data after HJD' = 980.0 have been re-reduced using image subtraction (Alard 
1999) which we found produces significantly lower errors and fails significantly less 
frequently than even the best DoPhot reductions previously reported by Albrow et 
al. (1999c). Details of this comparison will be given elsewhere. In addition, we have 
eliminated the SAAO data from the first night, HJD' = 975, because there was yet 
another CCD change on HJD' = 976.0 rendering the conditions on this night unique. 
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Finally, the fourteenth light curve is based on standard Johnson V band 
observations taken by the PLANET collaboration using the SAAO 1 m. These 
comprise 24 observations including 14 taken during the second caustic crossing and 
10 taken over the next 32 days. These data were used by Albrow et al. (1999a) 
to determine the V brightness of the source and so its V — I color, but have not 
previously been made available. 

In a preliminary fit to the data, we find that four data points are significant 
outliers. These are the SAAO point at HJD' = 979.6424, the MACHO B point 
at HJD' = 982.2061 the EROS R point at HJD' = 982.8427, and the MACHO 
B point at HJD' = 997.1607, with residuals of -5.4, -4.9, -4.1, and 3.9 a, 
respectively We eliminate these from future modeling. We renormalize the quoted 
errors from each light curve by a factor so as to force x 2 /dof (degree of freedom) 
to be unity for that curve. The factors are MACHO R (SoDoPhot): 1.12, MACHO 
B (SoDoPhot): 1.12, MACHO R (image subtraction): 1.26, MACHO B (image 
subtraction): 1.58, MACHO-CTIO R: 0.94, MACHO-CTIO B: 1.10, OGLE I: 
1.00, EROS R: 1.32, EROS B: 0.96, MPS R: 1.80, PLANET-SAAO (HJD' < 980) 
I: 1.04, PLANET-SAAO (HJD' > 980) I: 0.97, PLANET- Yale-CTIO I: 0.97, 
PLANET-CTIO I: 0.90, and PLANET-SAAO V: 2.21. The PLANET-Canopus 
I was not renormalized because there is only one point. The precise value of the 
renormalization factors depends slightly on which solution is adopted, but we find 
that our results are not sensitive to these small changes. We bin the early MACHO 
R and B data in 20-day intervals (see Fig. 1, below). With this binning, there are a 
total of 1018 data points. 

3. Method 

To analyze these data, we follow and slightly extend the method of Albrow et 
al. (1999c). We first review this method and its motivation and then discuss its 
extension. 

Events where a non-rotating binary lens passes in front of a uniform finite source 
are described by (7 + In) parameters, where n is the number of observatories: three 
parameters correspond to the three geometrical parameters of a point-like single 
lens (Einstein time scale t^, impact parameter uq, and time of closest approach to), 
three other parameters characterize its binary nature (mass ratio q, separation d 
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in units of the Einstein radius, and angle a of the source trajectory relative to the 
binary axis), one, p*, describes the size of the source relative to the Einstein ring, 
and there are n parameters for the source flux, F Si i . . . F Sin , and n for the unlensed 
background light F^i . . . F^ n , that is one pair for each of the n observatories. For 
events where the source crosses a fold caustic, one may define several additonal 
useful parameters which can be derived from these 7 + 2n, including the position 
within the Einstein ring where the source center crosses the caustic, u cc , the time 
of the caustic crossing, t cc , the angle <f> of the source trajectory relative to caustic 
at the crossing, the half-duration of the crossing At, and the radius crossing time 
= At sin 0. Generally, it is not difficult to find a set of parameters that yield a 
satisfactory fit to the data. However, it is often unclear whether there exist other 
equally good or better fits. One would like to make a systematic search through 
parameter space but because of the size and complexity of the parameter space, a 
brute-force search is out of the question. 

Albrow et al. (1999c) showed that for events with a well-observed caustic 
crossing, it is possible to greatly reduce the space of allowed solutions thereby 
permitting a systematic search of the remaining parameter space. The method 
proceeds in three steps. First, the caustic crossing is fit to a 5-parameter function. 
Second, these parameters are used to constrain a coarse-grained but systematic search 
through parameter space for solutions that can accomodate the non-caustic-crossing 
data. Third, final solution(s) are found by x 2 minimization using the results from 
the coarse-grained search as initial guesses. 

In the first step, the light curve is fit to a 5-parameter curve of the form 
where 

2 r 1 fl-x 2 \ 1/2 
G (v) = - dx [ LJL ) ( X 

7T Jmax(jj,-1) \ X — T] J 

is the normalized light curve of a (second) caustic crossing with a uniform source, and 
is a step function. Here Q is related to the rise time of the caustic (defined more 
precisely below), F cc is the magnified flux from the source when it is immediately 
outside the caustic, and u is the slope of the light curve immediately outside the 
caustic. Using the PLANET data, Albrow et al. (1999c) found 

Q = (15.73±0.35) F| day, t cc = (982.62439±0.00087) day, At = (0.1760±0.0015) day, 

(3) 



^cc) J 



-v), 



(2) 
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F cc = (1.378 ± 0.096) F 20 , Cu = (0.02 ± 0.10) F 20 day" 1 , (4) 

where F 2 o is the flux from an / = 20 star. For the first step, we simply adopt the 
results summarized in equations ([|) and (|). 

In the second step, the search of the full parameter space is substantially 
narrowed by making use of these caustic-crossing parameters with the following 
relations between observed and theoretical quantities, 

F C c = A CC F S + Ffo, -fbase = F s -\- Ffo, (5) 

t v = u T t E \ csc 0|, Q = F s 2 t r , (6) 

and 

At = t E p*\ esc 0| . (7) 

Here F s is the source flux, is the background flux, Fb ase is the baseline flux, A cc 
is the total magnification of the three non-divergent images at the position of the 
caustic, and u r characterizes the square-root singularity of the caustic. That is, in 
the neighborhood of the singularity, the total magnification of the two divergent 
images is given by A^{u) = (A-uj_/-u r ) _1 / 2 , where Au± is the perpendicular distance 
from the position u to the caustic in units of the Einstein radius, 9 E . The Einstein 
crossing time is t E — #e//^, where fi is the relative source-lens proper motion, 
6>| = (4GM/c 2 )(D LS /D L D s ), M is the total mass of the binary, D L and D s are the 
distances to the lens and source, and Dls = D s — D L . Finally, t T is a parameter that 
characterizes the rise time of the caustic. 

In an ideal world, t cc , Q, F cc , and Fb ase would be measured exactly. In this case, 
the search could be reduced to four parameters (d,q,£,t E ). (Recall that the ninth 
parameter, At, does not enter into the fit to the non-caustic-crossing data.) Here d is 
the binary separation in units of the Einstein ring, q is the binary mass ratio, and I is 
the position along the caustic of the second caustic crossing. For each pair (d, q), one 
steps along the caustic and determines A cc and u r from the lens geometry. Equation 
(|) yields F s = (F cc - F hase )/ (A cc - 1) and F h = F hase - F s . Next one chooses a value 
of t E . Equation (|6]) then fixes 0: | sin0| = F^u T t E /Q. This completely determines 
the geometry and source trajectory (up to a two fold ambiguity in 0). 

In practice, Q, F cc , and -Fb ase all have significant measurement uncertainties, and 
as Albrow et al. (1999c) discuss, this implies that one must allow a fifth continuous 
free parameter, 0, although this need be considered only over a restricted range. 
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If, as in the present case, light curve measurements come from several telescopes in 
several bands, then one must allow additional parameters for the source flux and 
background flux for each. However, these can be obtained from a simple linear fit 
(see also below) and so do not add significant computation time. The search through 
this 5-parameter space can be considerably simplified if there is information about 
the time of first caustic crossing. Then, for each trial trajectory, one can first check 
if the last measured point before the first caustic indeed lies outside the caustic and 
if the first point after the first caustic indeed lies inside. If either of these conditions 
is not satisfied, the trajectory can be rejected without further investigation. 

For the second step of the method, we follow Albrow et al. (1999c) with the 
following exceptions. First, we include the non-caustic-crossing data from all the 
light curves (except the SAAO V band data which are too sparse to contain useful 
information for this purpose). 

Second, we restrict the first caustic crossing to lie between the MPS point 
before the caustic at HJD' = 970.0503 and the OGLE point after the caustic at 
HJD' = 970.9037. As noted in § 2, the OGLE point could in principle be on the 
rising side of the first caustic crossing, i.e., before the center of the source crosses 
the caustic. However, a MACHO-CTIO R band data point taken approximately 1.2 
hours after the OGLE point rules out this possibility. This point was not reported 
by Alcock et al. (1999a) due to an oversight but was reported by Rhie et al. (1999). 
When the relative normalizations of the different light curves are properly set, the 
MACHO-CTIO point lies (20 ± 10)% below the OGLE point. If the OGLE point 
were before the first caustic, we would expect the light curve to be rising extremely 
rapidly, by of order a factor 2 in an hour, just as it is falling very rapidly at the end 
of the second caustic crossing (Afonso et al. 1998). Thus, the OGLE point certainly 
occurs on the falling side of the first caustic crossing. 

By restricting the first crossing to less than a day, we obtain a much more 
powerful constraint than the one used by Albrow et al. (1999c) who limited the first 
caustic crossing only to HJD' < 973.8. However, this change implies that smaller 
step sizes are required for and sin <fi so as to avoid missing the first caustic. We 
choose 2% increments for each compared to 5% used by Albrow et al. (1999c). Since 
the two caustics are separated by 12 days, there are guaranteed to be at least 3 time 
steps for which the first caustic crossing lies in the designated range. 

Finally, to avoid missing rotating wide binaries in the second step, we set the 
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model magnifications equal to unity (A = 1) for all points prior to HJD' = 810 (i.e., 
about 160 days prior to the first caustic crossing). The MACHO data are fairly flat 
during this period (Fig. [l]). In fact, while many of the binaries that we consider are 
at baseline during this early period, others, notably wide binaries, are not. They 
often show a "bump" (brightening then fading) several hundred days before the 
caustic crossing as the source approaches the companion star. Since this bump is not 
seen in the data, such binaries would seem to be ruled out. However, it is possible 
that the companion moved between the time that the source passed closest to the 
companion and the time when the source crossed the caustic (at which time the 
geometry of the event was primarily determined). If it moved sufficiently far during 
this interval, then the source would not have come close enough to the companion 
to cause a significant bump (see § 6.1). Thus, we include the early data (since it 
helps set the baseline) but do not allow it to rule out wide binaries until we have 
had a chance to examine the possibility that they might have avoided detection by 
rotating. 

From this step, we find two allowed regions in (d, q) space. One lies near 
(d, q) ~ (0.5,0.5), and the other near (d, q) ~ (3.5,0.4). Albrow et al. (1999a), 
Alcock et al. (1999a), Udalski et al. (1998b), and Rhie et al. (1999) all considered 
solutions in the general vicinity of the first region, but none considered solutions near 
the second. The two allowed regions lie in the lower right and upper left quadrants of 
the broad range of possible solutions shown in Figure 6 from Albrow et al. (1999c). 

For the third step, Albrow et al. (1999c) consider trial trajectories defined by 
seven parameters: the time t cc and the duration At of the caustic crossing, the 
Einstein time scale t^, the projected separation of the binary in units of the Einstein 
radius d, the binary mass ratio q = M 2 /M 1; the distance of closest approach (in 
units of the Einstein radius) of the source to the midpoint of the binary, u , and 
the angle a (0 < a < 2tt) between the binary-separation vector (M 2 to M\) and 
the proper motion of the source relative to the origin of the binary. (The center of 
the binary is taken to be on the right hand side of the moving source.) For each 
observation, the magnification is evaluated in one of two ways. If the source lies at 
least 3.5 source radii from the caustic, the magnification is simply the magnification 
of a point source. If it is closer, the finite size of the source is taken into account 
using the approximation 
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where u p is the position in the Einstein ring of the center of the source, and Au P) ±_ 
is the perpendicular distances from u p to the nearest caustic. If Au Pt ± > p*, then 
u q = u p . Otherwise, u g is taken to lie along the perpendicular to the caustic 
through Up and halfway from the caustic to the limb of the star that is inside the 
caustic. The perpendicular distance from u q to the nearest caustic is Au q ^, Al(u p ) 
is the magnification of the 3 non-divergent images at the position u p , A^iuq) is the 
magnification of the 2 divergent images at the position u q , and p* is the source size 
in units of the Einstein ring. See Figure 3 from Albrow et al. (1999c). The argument 
of Go is negative if the center of the source lies inside the caustic and positive if it 
lies outside. For each light curve i = 1, ... 14, we then use standard linear techniques 
to find the source flux F s ^ and background flux F^i that minimizes \\ > 

„ 2 _ [F S)i A fs (t k ) + F\i - F k f 

Xi = 2^ ~2 ' l y J 

k °k 

where F k and a k are the measured flux and error for the observation at time t k . (We 
follow Albrow et al. 1999c in constraining F Sti to be the same for the five PLANET 
light curves, % = 9, 10, 11, 12, 13, and in constraining F fe l0 = F b)13 .) 



3.1. Limb-Darkening Parameterization 



Equation (§) is valid in the approximation that there is no limb darkening. We 
model the surface brightness of the limb-darkened source by 



i - r 



3/ d 2 ^/ 2 



(10) 



where 9 is the angular position on the source star relative to its center, and V is the 
limb-darkening parameter. Note that with this formulation, there is no net flux in 
the T term, so F s remains the total flux. Convolving the T term with the square-root 
singularity of the caustic, we find the limb-darkened magnification is given by 



A(u p ) = A fs (u p ) + TA ld (u p ) A ld (u p ) = A° 2 (u q 



1/2 



Au 



p,_L 



(11) 



where H n (rj) = G n (rj) - G (t}), and 

(n + l)\ 



GM^TT- 1 ' 2 - 



1 _ T 2\n+l/2 

dx y , '— ,„ 8(1-7?). 



(n + 1/2)! J max ( v -i) [x-r]) 1 / 2 



(12) 
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Explicitly 

Gidv) = ~ E (3 - 2er/)(e - V ) 3 / 2 Q(e - rj), (13) 

£ = ±1 



where is a step function. To allow for limb darkening, we then modify equation 
(I): 

2 _ ^ [FsjA s fa) + ^jdjA fa) + ^Kj ~ Fk\ a\ 

k k 
The limb-darkening parameter for light curve % is then Y% = Fu^jF s ^. 

It is conventional to parameterize limb darkening by 



S(6) = 5(0) 



1 



e 2 \ 1/2 



(15) 



In this case, the flux associated with the limb-darkening term is not zero. Rather, 
it is a fraction (3/c — l)" 1 of the total flux. In a multi-parameter problem, the 
limb-darkening parameter then develops correlations with other parameters with 
which it has no physical connection. In our formulation, there is no net flux in the 
limb-darkening term, so the effect of limb darkening rapidly and explicitly vanishes 
far from the caustic crossing, 

Hi/2(V) - " Jo-(-^)~ 5/2 > fa « ( 16 ) 

Thus there are no spurious correlations. To make contact with the usual formulation, 
we note that 

(17) 



r + 2 



Limb darkening affects the magnification only if the source is transitting or is 
very close to the caustic. Thus, in principle it could affect the SAAO V and I curves 
(which both covered most of the caustic crossing), the Yale-CTIO curve (which has 
one point just before the end of the caustic crossing), the EROS B and R curves 
(which have 16 points each during the last 110 minutes of the caustic crossing), the 
MACHO CTIO R curve (which has 2 points just before the end of the crossing), 
and the MACHO B and R curves (which have points up to 1.7 At before the caustic 
crossing) . 

While the Yale-CTIO curve does not have enough coverage of the caustic 
crossing to make an independent estimate of the limb darkening, it is tied to the 
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SAAO photometry as discussed following equation (H) and more thoroughly in § 2 
of Albrow et al. (1999c). Thus, this one Yale-CTIO point can enter the fit for the 
SAAO / limb-darkening parameter. On the other hand, from equation ([16]), we see 
that limb darkening affects the MACHO B and R fluxes by less than a fractional 
amount YHi/ 2 {rj) /G (r)) ~ (3/160)r/~ 2 r <^ 0.3%, where we have assumed T <; 0.5 
and where we have made use of the limiting form Go(i]) ~ (— r/) -1 / 2 for i] <C —1. 
This compares to typical errors in MACHO photometry for these exposures of 2% 
to 3%. Hence we do not attempt to fit limb-darkening parameters to these two 
light curves. Thus, we fit for five independent limb-darkening parameters: SAAO 
V, EROS B, MACHO CTIO R, EROS R, and SAAO /, with corresponding central 
wavelengths of 0.55, 0.62, 0.64, 0.76, and 0.80 fim. 



4. Static-Binary Solutions 

We first search for solutions with static binaries. To do so, we will again set 
A = 1 for all points with HJD' < 810. In § 6, we will then investigate whether 
the solutions found in this way (or solutions near them) are in fact permitted 
when binary rotation is taken into account. We conduct the search on a grid with 
{Ad, Aq) = (0.02, 0.02) for the close-binary solution and (Ad, Aq) = (0.05, 0.04) for 
the wide-binary solution. 

We find two sets of static solutions. One is centered at (d, q) = (0.54,0.50). 
At the 3 a level (Ax 2 < 9), it extends from about (d, q) = (0.46,0.42) to about 
(d, q) = (0.60, 0.58), and is about half as wide in the orthogonal direction. The other 
solution is centered at (d, q) = (3.25, 0.24) and at the 3a level extends over the range 
d = 3.25ig;2o and q = 0.24tQ; 2 g. Dominik (1999b) has argued that there is a generic 
degeneracy in fitting light curves between a pair of close-binary and wide-binary 
solutions. The second (wide-binary) solution is formally favored at the 2 a level 
(A% 2 = 4), but we do not consider this to be a compelling reason to adopt it as the 
preferred solution. 

All the solutions near the close-binary minimum have similar parameters, as do 
all the solutions near the wide-binary minimum. For simplicity we quote the full 
set of parameters only at the minimum. These are shown in Tables 1 and 2. The 
division between the two tables is such that the parameters shown in Table 2 are 
derived from the linear fit described by equation (|T4"D and so have associated error 
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bars. The remaining parameters are shown in Table 1. Note that only the first 7 
parameters in Table 1 are independent. The five remaining parameters are derived 
from the fit. In particular, t* = At sin 0, and to = ^cc — ^(Mcc^cosa + Messina). 
We caution that the numbers of decimal places given for the parameters in Table 
1 convey a much higher precision than the statistical errors (which are in fact not 
even precisely known). The purpose of presenting many decimal places is to allow 
the reader to reproduce the solution. Because of strong correlations among the 
parameters, their values in a particular model must be known to high precision in 
order to avoid misdirecting the model into inappropriate regions of parameter space. 
The error bars in Table 2 reflect only the correlations within the linear fit described 
by equation ([14]) and not the correlations with the parameters in Table 1. Therefore 
these are actually lower limits on the errors. Note that we show the ratio F^/F s only 
for the light curves for which it is reasonably well determined (< 10%). 

The parameter that varies the most over the allowed set of the solutions is 
tE, which ranges from about 75 to about 125 days within the 3 a range of the 
close-binary solutions and from about 145 to 200 days within the 3 a range of the 
wide-binary solutions. From the standpoint of the proper-motion measurement, 
three parameter combinations are important, (V — I) s , I s , and t* = At sin 0. The 
first essentially does not vary at all, (V — I) s = 0.30 for all allowed solutions. When 
(V — I) s is fixed, the proper motion scales as \i oc lO -0,27 *^" 1 . The full (3cr) range of 
variation of this parameter combination and thus of \x is only about 25% over each 
the two classes of solutions. 

The limb- darkening coefficients given in Table 2 are shown in Figure |2|. The 
close-binary is shown by open circles and the wide-binary is shown by filled circles. 
The horizontal error bars show the full-width at half maximum (FWHM) of the 
filters, while the vertical error bars denote the statistical errors. We emphasize 
again, however, that these include only the errors from the linear fit generated by 
equation (|T4| ) and not those that arise from correlations with the other parameters. 
If we suppress limb darkening and force a fit to a uniform disk, then in both cases, 
X 2 increases by about A% 2 = 38 for 5 degrees of freedom. This is far less than the 
A% 2 = 106 that would be predicted based on a naive interpretation of the errror bars 
shown in Figure 0, and this difference arises exactly from the fact that these error 
bars do not account for the correlations with other parameters. Nevertheless, the 
full fit reveals that limb darkening has been detected with high significance (formally 
99.999%). 



-17- 



Unfortunately, to the best of our knowledge, no theorists have ever calculated 
limb-darkening profiles of metal-poor A stars. Since limb darkening has clearly 
been detected in one such star, perhaps some theorist will now undertake such a 
calculation. For the close-binary solution, the limb-darkening coefficients fall from 

0. 45 ± 0.11 for V to 0.17 ± 0.04 for I. The wide-binary solution is similar. The 
one exception is MACHO-CTIO R, but its error bars are too large to make any 
definite statement because its limb-darkening parameter was derived from only two 
measurements. 

Tables 1 and 2 also show a third solution, one for a rotating wide binary. This 
solution is derived in § 6.1, below. As is clear from Figure [I], the static wide-binary 
solution is not viable: it is only an intermediate step on the way to finding a viable 
rotating wide-binary solution. Hence, in Tables 1 and 2, the close-binary and 
rotating wide-binary solutions are labelled "viable" while the static wide-binary 
solution is labelled "not viable". However, until we introduce rotation in § 6, all 
references to the wide-binary solution will be to the static version. 

Figures ^| and ^ show the model light curves together with all the available data 
for the close-binary and wide-binary solutions respectively. Because the data are in 
different passbands, we cannot compare the predicted flux with the observed flux 
as we could if the data were in a single passband. We therefore deblend our data, 

1. e. we plot 2.5 log[(F — F h )/F s ] (points) and compare this to 2.5 log (magnification) 
(solid curve), where F s and F h are the fit values of the source and background flux. 
The points are binned primarily in one-day intervals. However, the points before 
HJD' = 950 are binned in 10-day intervals and the points near the caustic crossings 
are binned in 0.1-day intervals. Data from different observatories are combined 
together. Figures [5] and |] show close-ups of the two model fits in the neighborhood 
of the second caustic crossing binned in 0.01-day intervals. 

The two fits appear to be equally good to the eye. This is illustrated in Figure |7| 
which shows the fractional difference in the predicted fluxes between the two models 
for each of the 14 light curves analyzed in this paper. The fundamental physical 
reason for this degeneracy is shown in Figure |8| where the caustic structures for the 
two solutions are superposed. These caustic structures are very similar. 
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5. Proper Motion 

The proper motion is given by /i = To obtain the proper motion one must 

therefore estimate the angular source size 9*, which can be calculated if one knows 
the dereddened color and magnitude of the source. Among all the light curves, there 
are photometric calibrations for only five: PLANET V (SAAO only) and PLANET 
/ (SAAO and Yale-CTIO) (Albrow et al. 1999a), MACHO B and R (Alcock et 
al. 1999b), and OGLE /. As we describe below, the calibration of PLANET / 
is tied to the OGLE calibration. We find that the F s values for these two light 
curves are consistent at the 1 a level. However, the errors for the OGLE I F s are an 
order of magnitude larger than for PLANET / (because there are many fewer data 
points), so the OGLE / F s does not yield significant additional information about 
the flux of the source. For the close-binary solution, we have V s = 22.36 ± 0.01, 
I s = 22.06 ± 0.00, (V - I) s = 0.30 ± 0.01, from PLANET and V s = 22.67 ± 0.01, 
R s = 22.58 ± 0.01, (V - R) s = 0.09 ± 0.01 from MACHO. For the wide-binary 
solution, we have V s = 22.17 ± 0.01, I s = 21.87 ± 0.00, (V - I) s = 0.30 ± 0.01, from 
PLANET and V s = 22.48 ± 0.01, R s = 22.39 ± 0.01, (V - R) s = 0.09 ± 0.01 from 
MACHO. 

In addition to these errors reported by the fit, Albrow et al. (1999a) estimate 
that their calibration error is 0.02 mag for the PLANET color and Alcock et 
al. (1999b) estimate that their calibration error is 0.04 mag for the MACHO 
color and 0.10 for the magnitudes. Two points are clear from this summary. 
First, the ratios of fluxes are essentially identical for the two models. Second, 
the MACHO and PLANET colors are mildly inconsistent and the MACHO and 
PLANET V magnitudes are inconsistent at the 3<r level. We believe that the 
PLANET calibration is substantially more reliable than the MACHO calibration 
since PLANET calibrated their data using secondary standards in the field that were 
in turn measured in the standard way by OGLE (Udalski et al. 1998a), i.e. from 
primary standards on photometric nights. On the other hand, although MACHO 
applies essentially the same procedure for their calibration of their LMC fields, for 
the SMC they simply adopt the mean zero points derived for the LMC at similar 
airmass (Alcock et al. 1999b). We therefore adopt the PLANET calibration. 

Following Albrow et al. (1999a), we adopt a total extinction of Ay = 0.22 ±0.1. 
The final results do not depend strongly on the extinction (see below). The flux is 
given by F = 9^S, where S is the mean surface brightness of the source. We will 
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assume that this surface brightness is a function only of the (V — I)o color and not 
any other properties of the star. (We know, for example, that the star is a dwarf 
rather than a giant.) We can then write 

6* = 79 nas KT - 2 ^- 22 ) ( — — ) (18) 

\D(y_/) 0= o.2l/ 

where we have evaluated the normalization using Green, Demarque, & King (1987), 
specifically their Y = 0.2, Z = 0.001, Age = 1 Gyr table. We therefore obtain 
estimates of 82 nas and 89 nas for the angular size of the source in the close-binary 
and wide-binary solutions respectively. As described by Albrow et al. (1999a), this 
estimate has a 3% error for uncertainty in the extinction Ay (Albrow et al. 1999a) 
and a 5% error for uncertainty in the theoretical model (M. Pinsonneault 1998, 
private communication), for a total uncertainty of 6%. 

Hence in the two models the proper motions are 

H = 1.30 ± 0.08 km s" 1 kpc" 1 (close binary) (19) 

and 

fi = 1.76 ± 0.11 km s _1 kpc -1 (wide binary) (20) 

The errors in these equations reflect only the uncertainties in the extinction and the 
stellar models, and they do not include uncertainties in the parameter fits. Recall 
from § 4, however, that even at the 3 a level, the range of allowed values of the 
parameter combination lO -0 " 2- ^" 1 oc \i is very restricted. 

The values in equation (|T9| ) and ( p0[ ) clearly put the lens in the SMC rather 
than the Galactic halo. For comparison note that Albrow et al. (1999c) found some 
solutions that were moving much faster and hence would not be explainable as SMC 
events. These additional solutions are ruled out by combining all the available data. 

5.1. Binary Physical Characteristics 

Since the binary is known to be in the SMC, we can use the proper-motion 
measurements to obtain estimates of the binary physical projected separation, 



r p = <i/it E D s ~ djjLt-^ x 60kpc. 



(21) 



-20- 



This yields r p = 2.40 AU and r p = 32.4 AU for the close-binary and wide-binary 
solutions respectively. Note that for circular face-on orbits, r p = a, the semi-major 
axis, while all orbits satisfy a > r p /2. 

These projected separations can be used to place limits on the motion of the 
binary. For example, for the close-binary solution, the blended background flux in 
the MACHO B band light curve (which has the best-determined blended flux and is 
also the only well-determined measurement in the blue) is approximately 50% larger 
than the source flux, so the larger of the two lens stars cannot be more than about 
2.5 Mq. For the wide-binary solution, the blended and source fluxes are about equal, 
so the larger star cannot be more than about 2M . Thus the total mass of the 
binary in both cases is limited to M ^ 3M . If we momentarily assume a face-on 
circular orbit, then from Kepler's Third Law, the period is constrained to P > 2yr 
and P > HOyr for the two solutions. For a face-on eccentric orbit at apastron, 
the periods could actually be shorter by 8 1 / 2 , but what actually concerns us is not 
the length of the period but the relative motion of the binary lenses over times 
that are very short compared to the period. For a circular orbit, this instantaneous 
angular speed is u C i TC = 27r/P circ , but the maximum instantaneous angular speed 
occurs for a face-on eccentric orbit where the caustic crossing occurs near periastron: 
w m£H = 2 1 / 2 u; circ . We must therefore consider binary motions up to this level. 

6. Rotating Binaries 

Although binaries are not static, only a few attempts have been made to fit 
microlensing light curves to dynamic binaries (Dominik 1998). In principle, it is 
possible to measure six orbital parameters of a binary from sufficiently precise 
observations. These are the same six that can be measured from proper-motion 
measurements of visual binaries except that the angular semi-major axis is measured 
relative to #e (rather than absolutely) and the line of nodes is measured relative 
to the direction of the source (rather than celestial coordinates). In practice, it is 
extremely difficult to measure anything other than the (two-dimensional) projected 
relative velocity of the components in units of $e- In fact, no binary-motion 
information of any type has ever been extracted from a microlensing event. We 
therefore restrict consideration to the simplest form of such motion, uniform circular 
motion in the plane of the sky. This leaves the geometry of the lens fixed and 
permits only rotation of this geometry. If we allowed more general two-dimensional 
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motion, the geometry of the lens would change as the projected positions of the two 
components moved closer together or further apart. We will explicitly ignore this 
type of change in the binary configuration. 

6.1. Wide-Binary Solution 

As we discussed in § 3, we forced the magnification at early times to A = 1 
when fitting the static binary solutions because the MACHO light curve is observed 
to be flat at these times. Had we not done so, the wide binary solution would have 
been ruled out at the 18 a level (A% 2 = 342). In Figure [I] we show the early light 
curve for the best-fit static-binary solution together with the MACHO data. The 
model is clearly ruled out by the data. In fact, we find that even if we allow this 
binary to rotate with a period of P = 75 yr, i.e., the minimum permitted by the 
argument of § 5.1, the model is still ruled out at the 8 a level (Ax 2 = 58). However, 
there are satisfactory rotating binaries in the neighborhood of the best-fit static 
solution. In Tables 1 and 2 we give the parameters for one of these rotating solutions 
with (d, q) = (3.65, 0.36) and P = 75 yr, and in Figure |^ we show a diagram of the 
caustic structure for this rotating solution together with the corresponding static 
solution for the same (d, q). Since the two lenses are separated by much more than 
an Einstein ring, the magnification structure is for the most part a superposition 
of the magnification of two isolated lenses. Hence, it is clear from the diagram why 
the static model is excluded: the source passes within ~ 0.4 binary-mass Einstein 
radii of the larger lens, which is about 0.55 Einstein radii scaled to the mass of 
this lens. Thus, the magnification is about 2. Even though the event is heavily 
blended (F^/F s ~ 2) and the errors in MACHO photometry are relatively large at 
these early times, this magnification would easily be seen in the data. However, if 
the binary were rotating clockwise in the plane of the sky, then two years before 
the caustic crossing at the time of closest approach, the source would have been 
about twice as far from the heavier lens, thus reducing (A — 1) by a factor 3.5. We 
show the light curve resulting from this rotating-binary solution in Figure [I]. It 
is barely distinguishable from the baseline. Numerically we find that the 75-year 
period binary increases y 2 by less than 1 unit relative to the artificial case used in 
the initial simulations of A = 1 for HJD' < 810. We find that the y 2 of this rotating 
wide binary is only 5 higher than the \ 2 °f the best-fit close binary. We conclude 
that the data are consistent with a wide-binary solution. The upper panel of Figure 
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H is a close up of the trajectories with and without rotation. The two trajectories 
are essentially identical in the region around the caustic crossing. 

Because the allowed rotating binaries tend to be on one side of best-fit static 
wide-binary solution (higher d and higher q), they tend to have systematically lower 
proper motions than that of the static solution given in Tables 1 and 2. For the 
rotating wide binary shown in Figure [| (and indeed for its static analog), we find 

ji = 1.48 ± 0.09 km s -1 kpc -1 (rotating wide binary). (22) 

6.2. Close-Binary Solution 

For the close-binary solution the rotation periods can be much shorter, so that 
in contrast to the situation illustrated in Figure |^ for the wide binary, the rotating 
and non-rotating trajectories are not the same within the caustic region. Hence, 
rotating solutions require substantially different geometries. For example, for P = 10 
yr (counter-clockwise), we find a best fit at (d, q, t^j = (0.56,0.44,95.5 days), and 
for P = 10 yr (clockwise), we find (d,q,t-E) = (0.54,0.58,95.9 days). These solutions 
are about 1 a worse fit than the non-rotating solutions. Their proper motions 
are respectively 4% lower and 7% higher, and the limb-darkening parameters are 
very similar to the non-rotating case. Thus, while both rotating and non-rotating 
solutions are compatible with the data and while allowing rotation increases the 
uncertainties in the binary parameters, these various solutions have very similar 
implications for the nature of the source and lens. 

7. Comparison With Previous Solutions 

Six papers have made estimates of some or all of the parameters of MACHO 
98-SMC-l based on subsets of the data presented here (Afonso et al. 1998; Albrow 
et al. 1999a,c; Alcock et al. 1999a; Udalski et al. 1998b; Rhie et al. 1999). We 
now analyze the relationship of the results presented in this work to these earlier 
efforts. We will focus attention on whether the various previous solutions and partial 
solutions are consistent with one another and with the present results, and we will 
attempt to resolve any inconsistencies. 

Rhie et al. (1999) analyzed almost all of the non-PLANET data presented here, 
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and Albrow et al. (1999c) analyzed almost all of the PLANET data. Hence, our 
present analysis is essentially based on the union of these two disjoint data sets. 

Albrow et al. (1999a) published two solutions, now known as PLANET Model 
I and PLANET Model II. However, Albrow et al. (1999c) subsequently showed 
that PLANET Model I actually sits in a single extremely broad, virtually flat, x 2 
minimum which connects all of the close-binary solutions that they found. PLANET 
Model I is essentially the same as model #26 from Albrow et al. (1999c). Because 
of their excellent coverage of the second caustic crossing, Albrow et al. (1999c) were 
able to measure t cc very precisely and At fairly precisely. Those measurements are 
confirmed by the solutions presented here. On the other hand, they showed that the 
broad degeneracy in their overall solution could be traced to their lack of coverage 
of the early light curve (see their Figs. 7, 8, and 9). One would expect as more 
data are added to the data available to Albrow et al. (1999c), that the two broad 
minima shown in their Figure 6 would contract and possibly break up into several 
discrete local minima. The close-binary solution presented here is very similar to an 
interpolation between the neighboring grid points of their models #27 and #31. 

Afonso et al. (1998) measured the parameter combination t cc + At = 
982.8039 ± 0.0010 (after correction of a transcription error in the original paper) 
based on EROS coverage of the end of the light curve. This differs by only 2 
minutes from the values shown in Table 1. Alcock et al. (1999a) modeled the 
event by combining their own MACHO/GMAN data with the EROS data. The 
MACHO/GMAN model was refined by Rhie et al. (1999) after the time of the 
first caustic was pinned down by their own MPS data together with the OGLE 
data (Udalski et al. 1998b). We now investigate the consistency of the MPS model 
with the models based solely on the PLANET data on the one hand, and with the 
close-binary (CB) model presented here on the other. 

The first question to ask is: are the MPS and CB models in discrete local 
minima, or are they two different points in the same minimum? They are located 
at (d,q,t E ) = (0.646,0.518,70.5) and (0.54,0.50,99.0) respectively. To answer this 
question, we find solutions based on all the data, but subject to the constraint 
of fixed {d,q). We evaluate these solutions on a grid of (Ad, Aq) = (0.02,0.02) 
in the neighborhood of the CB model at (d, q) = (0.54,0.50). We find that x 2 
varies smoothly over this grid of solutions. The solution at (d, q) = (0.64, 0.52) is 
extremely similar to the MPS solution, and x 2 rises monotonically between the CB 
and MPS-like solutions (Ax 2 = 32). Hence these two solutions are in the same 
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minimum and are not discrete minima. 

Since the MPS solution has higher x 2 based on all the data and is in the same 
minimum as the CB solution, did MPS therefore find a false minimum? To address 
this we evaluate x 2 f° r the two solutions based on the data available to MPS (i.e. 
excluding the PLANET data and the EROS data from other than the night of the 
caustic crossing and using the SoDoPhot reductions of the MACHO data rather 
than image subtraction.) We then find that the MPS solution is favored over CB by 
Ax 2 = 10. That is, the two solutions differ because they are based on different data 
sets rather than because of different modeling procedures. 

Finally, we investigate a conflict between the MPS and PLANET solutions 
which was previously identified by Rhie et al. (1999). They noted that their value 
of = 982.683 or t cc = 982.694 (depending on their model of limb darkening) is 
later than the PLANET value t cc = 982.62439 ± 0.00087 (Albrow et al. 1999c and 
confirmed here). They did not quote error bars on their own value which is based 
on modeling the interpolation between MACHO data cutting off 0.3 days before the 
crossing and EROS data beginning 0.1 days after it. However, by evaluating x 2 for a 
series of models with the parameters d, q, a, uo, £e, and t cc + At fixed, but t cc varying, 
we find that the error in the MPS value for t cc is approximately 0.009 days. Thus, 
the difference between the MPS and PLANET values is a 6 a discrepancy if due to 
an MPS problem and a 68 ex discrepancy if due to a PLANET problem. Clearly this 
difference is not the product of a statistical fluctuation. 

To determine the origin of this conflict, we fit all the data but use the original 
SoDoPhot reductions (used by MPS) in place of the image-subtraction reductions 
(used here) for the MACHO data. We find that the MACHO data points lie 
systematically above the model at the beginning of night before the caustic crossing 
(HJD' ~ 982.1) and systematically below the model at the end of the night 
(HJD' ~ 982.3). There is no such systematic trend in the MACHO data points when 
reduced by image subtraction. We infer that this trend may have been responsible 
for the late t cc in the MPS model. We test this hypothesis by redoing the fits 
based on an MPS-like data set but substituting image-subtraction reductions for 
SoDoPhot. We then find that the (d, q) = (0.54, 0.50) solution is favored over the 
(d, q) = (0.64, 0.52) solution by Ax 2 = 5. Moreover the t cc in the best-fit model now 
differs from the PLANET value by only 0.02 days which is only a 2 a discrepancy. 

In brief, the EROS measurement of t cc + At and the PLANET measurement 
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of t cc have been confirmed with high precision. The original MACHO model when 
refined by MPS based on the MPS+OGLE determination of the first caustic crossing 
holds up very well. It lies close to the CB solution based on all the data. In 
hindsight, image subtraction would have yielded even more precise refinements of 
this model. Finally, one sub-region of the broad class of wide-binary solutions found 
by PLANET (Albrow et al. 1999c) survives the inclusion of the non-PLANET data. 
We pat our collective selves on the back for a job well done. 



8. Conclusions 

We have combined the data on MACHO 98-SMC-l from five collaborations 
to produce one of the best sampled microlensing light curves ever published. We 
confirm earlier claims that the relative source-lens proper motion is low, so the lens 
must be in the SMC. However, there is a twist: despite the fact that our combined 
data set is enormously superior to any of the individual data sets, there are two 
very distinct solutions that are compatible with all the data, a close-binary and a 
wide-binary solution. Fortunately, both have very similar proper motions so there 
is no significant ambiguity in this parameter. We find a relative proper motion of 
/i ~ 1.30 kms -1 kpc~ x or /i ~ 1.48 kms -1 kpc -1 . 

We have measured the limb-darkening parameter in five different bands with 
centers at 0.80, 0.76, 64, 0.62, and 0.55 /im. If our results are expressed in terms of 
the standard limb-darkening parameter c, the respective values for the close-binary 
solution are 0.23 ± 0.05, 0.24 ± 0.05, 0.06 ± 0.34, 0.42 ± 0.05, and 0.55 ±0.11. All 
other solutions have limb-darkening parameters that are close to these. 
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Fig. 1— MACHO B and R data for MACHO 98-SMC-l binned in 20-day intervals 
for the period before HJD' = 810. The bold lines are the values for the baseline 
flux from the close-binary solution, (d, q) = (0.54,0.50). The event shows no 
significant deviation from baseline during this early period. The solid lines are from 
the best fit for the non-rotating wide binary (d, q) = (3.25, 0.24) which is clearly 
ruled out by the data. However, a wide binary with a 75 year period and a nearby 
(d, q) = (3.65, 0.36) is permitted. Its early light curve is shown as a dashed curve that 
is barely distinguishable from a flat-baseline curve. 



Fig. 2.— Limb- darkening parameters T (as defined in eq. |TUp derived from the 
close-binary solution (d, q) = (0.54,0.50) (open circles) and the wide-binary solution 
(d, q) = (3.25, 0.24) (filled circles) to MACHO 98-SMC-l for 5 passbands. From left 
to right: PLANET-SAAO V, EROS 5, MACHO-CTIO R, EROS R, and PLANET- 
SAAO /. The horizontal error bars represent the full-width at half-maximum of the 
filters, and the vertical error bars are statistical. The limb- darkening parameters 
for the two solution are very similar because the measurement of limb darkening 
depends primarily on the caustic crossing and not on the global characteristics of 
the light curve. To avoid clutter, the error bars for the wide-binary solution are not 
shown, but they are almost identical to the error bars for the close-binary solution. 



Fig. 3. — Predicted versus "observed" deblended magnification for the close-binary 
model (d, q) = (0.54,0.50). The deblended magnification is A = (F — F h )/F s where 
F is the observed flux, and F s and F h are the fit source and background fluxes in 
the model. Data are binned, mostly in 1-day bins. However, for HJD' < 950 there 
are 10-day bins, and in the immediate neighborhood of the caustics there are 0.1-day 
bins. Data from all 14 light curves from the 5 collaborations are averaged together 
whenever they lie sufficiently close to fit in the same bin. 



Fig. 4. — Predicted versus "observed" deblended magnification for the wide-binary 
model (d,q) = (3.25,0.24). Similar to Fig. 0. 



Fig. 5. — Predicted versus "observed" deblended magnification for the close-binary 
model (d, q) = (0.54, 0.50) showing the vicinity of the second caustic crossing. Same 
as Fig. |3| except that bins are 0.01 days. 
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Fig. 6. — Predicted versus "observed" deblended magnification for the wide-binary 
model (d, q) = (3.25, 0.24) showing the vicinity of the second caustic crossing. Same 
as Fig. £| except that bins are 0.01 days. 



Fig. 7. — Fractional differences between fluxes predicted by the close-binary solution, 
(d,q) = (0.54,0.50), and the wide-binary solution, (d,q) = (3.25,0.24), for the 14 
different light curves (solid lines). 



Fig. 8. — Caustic structures for the close-binary (bold curve) and wide-binary (solid 
curve) solutions. Each has been rescaled according to the Einstein crossing time of 
the solution. The caustics have been rotated so that the source trajectories (straight 
solid line) overlap. Time is shown in days from the second caustic crossing, so source 
motion is to the right. 



Fig. 9. — Wide-binary trajectory with (d, q) = (3.65, 0.36) for static case and for 
binary with P = 75 yr period. The upper panel shows a close-up of the caustic 
together with the two source trajectories which are barely distinguishable. The light 
curve in this region is therefore independent of rotation and the structure of the 
caustic fixes the local trajectory. The lower panel shows the full caustic structure. 
The two caustics are separated by about 3.65 Einstein radii, or about 2 years. This 
interval is sufficient to allow the source-companion closest approach to grow by by a 
factor ~ 2 relative to the static case. This in turn reduces (A — 1) a factor of ~ 3.5. 
The static solution is ruled out by the data but the rotating solution is permitted 
(Fig. [I]). The source is moving to the left. 
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Parameter 


Close Binary 
Viable 


Wide Binary 
(Static) 
Not viable 


Wide Binary 
(Rotating) 
Viable 


d 


0.54 


3.25 


3.65 


Q 


0.50 


0.24 


0.36 


a 


350.575° 


173.344° 


172.922° 


u 


0.045033 


0.172928 


0.211577 


*e 


98.956 


164.728 


198.323 


tec 


982.62408 


982.62389 


982.62414 


At 


0.17836 


0.17889 


0.17880 




36.9° 


28.8° 


29.1° 




0.1071 


0.0862 


0.0869 




0.146824 


-1.412225 


-1.652962 


^cc,y 


0.021277 


-0.009302 


-0.007959 


to 


968.63593 


751.73636 


657.49642 


x 2 


981.1 


976.8 


986.4 



TABLE 2 



Parameter 


Close Binary 


Wide Binary 


Wide Binary 






(Static) 


(Rotating) 




Viable 


Not viable 


Viable 


T (SAAO /) 


0.17 ±0.04 


0.15 ±0.04 


0.15 ±0.04 


T (EROS R) 


0.17 ±0.04 


0.15 ±0.04 


0.15 ±0.04 


T (CTIO R) 


0.04 ±0.23 


0.01 ±0.26 


0.04 ±0.25 


T (EROS B) 


0.33 ±0.04 


0.31 ±0.04 


0.31 ±0.04 


T (SAAO V) 


0.45 ±0.11 


0.41 ±0.11 


0.40 ±0.11 


v s 


22.356 ±0.007 


22.171 ±0.007 


22.527 ±0.007 


Is 


22.058 ± 0.003 


21.875 ± 0.003 


22.231 ±0.003 


(V - i) s 


0.298 ±0.007 


0.297 ±0.007 


0.296 ±0.007 


F h /F s (MACHO R) 


1.98 ±0.04 


1.51 ±0.03 


2.40 ± 0.04 


F h /F s (MACHO B) 


1.57 ±0.03 


1.16 ±0.03 


1.92 ±0.03 


Fh/Fs (CTIO R) 


0.63 ±0.07 


0.65 ±0.06 


1.01 ±0.08 


F h /F s (SAAO /) 


1.48 ±0.06 


1.14 ±0.05 


1.74 ±0.07 


\i (km/s/kpc) 


1.30 ±0.08 


1.76 ±0.11 


1.48 ±0.09 
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